clc;
clear;
load dataAnalysisQ_G1000k.mat;
load MacroData_Q; % variable name is data


% 1 order difference
order1Diff = diff (data);
x1 = order1Diff (4:size(order1Diff,1), 3: size(order1Diff,2));
y1 = data(5:size(data,1), 2);
[r,beta]=ols_multi( y1, x1 ,x1 );
out=fsra(y1 ,x1);

% detrend
detrend_x = detrend(data);
x2 = detrend_x (4:size(detrend_x,1) , 3: size(detrend_x, 2));
y2 = y1;

% hp filter
[T,C] = hpfilter(data);


